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Abstract 

We consider some aspects of counting BPS operators which are annihilated by two 
supercharges, in superconformal field theories. For non-zero coupling, the correspond- 
ing multi-variable partition functions can be written in terms of generating functions 
for vector partitions or their weighted generalisations. We derive asymptotics for the 
density of states for a wide class of such multi- variable partition functions. We also 
point out a particular factorisation property of the finite ./V partition functions. Finally, 
we discuss the concept of a limit curve arising from the large ./V partition functions, 
which is related to the notion of a "typical state" and discuss some implications for 
the holographic duals. 
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1 Introduction 

Understanding the microscopic origin of black hole entropy is an important problem with the 
potential to shed light on Quantum Gravity. String theory has made remarkable progress 
in this area; for supersymmetric asymptotically flat black holes in four and five dimensions 
we know how to count the microscopic states in terms of D-brane degrees of freedom [1, 2]. 
However, an important factor which simplifies the counting in all known examples is the 
existence of an effective string picture - the microscopic theory in the regime of interest can 
always be modeled by a 1 + 1 dimensional CFT with an effective central charge [3]. The 
degeneracy of states is then simply obtained by applying the Cardy formula. With the 
recent advances in constructing supersymmetric black holes, we now have examples where 
this effective string model is insufficient; for example, the supersymmetric AdSs black holes 
[4, 5, 6, 7, 8] 1 are not described by an effective string model. To obtain a microscopic 
derivation of the degrees of freedom for such black holes, one needs to directly count the 
states (or operators) of the dual M = 4 super Yang-Mills (SYM) theory. 

The counting of operators in M = 4 SYM was initiated in [11, 12] and later extended to 
supersymmetric operators in [13]. More generally, one can consider the problem of counting 
chiral operators in Af = 1 theories (these are analogous to counting §-BPS states in J\f = 4 
SYM) . There is by now a large literature on counting such operators from various perspectives 
such as Polya counting, giant gravitons, dual giants, the 'plethystic program', character 
formulae, linear sigma models, etc., - see [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 
27, 28, 29]. 2 In M = 1 SCFTs one always has an R-symmetry, and eigen-operators of this 
symmetry must have discrete quantum numbers. The conformal dimension of such operators 
must then satisfy A > 3R/2, with equality if and only if the operator is BPS (i.e., annihilated 
by a supercharge). In addition, such theories will typically possess global symmetries and 
therefore operators will be labeled by quantum numbers corresponding to these (e.g., M = 4 
actually has an SU(A) symmetry which contains U(l) 3 ). Thus generically, BPS operators 
in such theories will be labeled by a set of k integers rij, with the conformal dimension A 
determined by some combination of them (corresponding to the R-symmetry). It is then 
found that in the large N limit one typically encounters partition functions of the form 

3*G9i, n \ — 

C" n 1 — e ^*=i PiUi 

71; >0 

The (3i are some set of chemical potentials conjugate to the quantum numbers in ques- 
tion rij. Given such a partition function one would like to extract the density of states, 

^ee also [9, 10] for some progress towards determining the most general such black hole via a classification 

of their near-horizon geometries. 

2 One can also consider BPS states in the D = 3 M2-brane SCFT and the D = 6 M5-brane SCFT [30, 31]; 

our discussion incorporates these cases as well. 
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d k (ni, n 2 , • • • n k ) defined by 

3*(/3i,/V--&) = J2 d k (n 1 ,n 2 ,---n k )e-ZUPim ( L2 ) 

ni,n 2 ,...,n k >0 

at least in the limit of large rij. 

In the k = 1 case the partition function (1.1) is the generating function for the number of 
partitions of an integer and corresponds to the classic result due to Hardy and Ramanujan: 3 

d\{n) — p{n) ~ pexp I ir\l —n J (1.3) 




as 4 n — > oo. Similarly, for k > 1 it turns out that (1.1) is the generating function for 
the number of partitions of the "vector" (ni, n 2 , ■ ■■, n k ) (we will define this precisely below). 
Curiously, it appears that asymptotics for the number of partitions of a vector are hard to 
come by in the literature. In [24] leading asymptotics for log d k (rii, n 2 , ■■.rik) for large rii (of 
the same order) were determined. 5 One of the main purposes of this paper is to improve 6 
upon this and derive the leading asymptotics for dk(ni,ri2, ■■■rik) when all the rij are large 
(and of the same order). Thus our results will be direct analogues of (1.3) for the number 
of vector partitions. 

A technical objection to determining such asymptotics would be that such multivariable 
partition functions do not seem to have any obvious modular properties, analogous to the 
single variable case (although see [33] for modular properties of free massless scalars in 
higher dimensions). Indeed the k = 1 case due to Hardy and Ramanujan can be deduced 
from the modular properties of the Dedekind function. However, one can also obtain (1.3) 
using a method due to Meinardus [34, 35], which does not employ any modular properties, 
and indeed we will show that the same method can be used to analyze these multivariable 
partition functions. 

For k > 1 the special case = (3 of the above partition function simplifies to 

oo 1 



n=l 



which is a generating function for c k {n) = ^2 m+n2+ _„ +nk=n d k (ni,n 2 , ...,n k ). At first sight it 
seems that the asymptotics in this special case can be deduced from the theorem of Meinardus 



3 Hardy and Ramanujan in fact obtained a much more precise formula which was improved upon later by 

Rademacher into an exact formula for p(n) . 

4 In this paper we use the mathematician's meaning of ~; that is, f(n) ~ g{n) as n — > oo is equivalent to 

5 See also [32] for the leading exponential behaviour of the k — 2 partition function encountered in the 
plane wave matrix model. 

6 Observe that if log<i(n) ~ g(n) it does not follow that d{n) ~ e 9 ^ n \ The converse, though, is of course 
true. 
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[34], a direct generalisation of (1.3) which applies to generating functions of the form (1.4) 
for some class of sequences a n . The class of sequences a n to which the theorem applies is 
encoded in the behaviour of the associated Dirichlet series 



In particular it is restricted to situations where the Dirichlet series D(s) has one simple pole. 
It is straightforward to show that the Dirichlet series defined by the a n above (1.4) actually 
has k simple poles. It should be noted that if one simply requires the leading asymptotic for 
logCfc(n), then one can use a trick employed in [17, 24, 36]. This involves writing (1.4) as 
a product of k generating functions whose associated Dirichlet series each have one simple 
pole, and observing that the leading asymptotic for logCfc(n) comes from the Dirichlet series 
with the rightmost pole (this is equivalent to the rightmost pole of (1.5)). However, this 
trick does not seem to generalise to determining the asymptotic behaviour of c(n). For this 
we return to the proof of Meinardus' theorem and consider associated Dirichlet series with 
multiple poles. This is not a fundamental obstacle and can be straightforwardly accounted 
for; therefore we obtain the leading asymptotic behaviour of Cfc(n) as n — > oo. 

In addition we discuss aspects of the finite N partition function. This is based on the fact 
that generating functions for the finite N partition function have been written down. These 
partition functions usually take the form of a grand canonical partition function (where N is 
the particle number) for bosons in a /c-dimensional harmonic oscillator potential. Combina- 
torially they can be interpreted in terms of restricted vector partitions, i.e., into partitions 
with no more than N parts. By exploiting the combinatorial properties of the partition 
functions we easily write down explicit expressions for the finite N partition functions and 
prove an interesting factorization property. 

We also discuss how one can use the data in the partition sum to talk about typical states 
in the ensemble. The latter while being an interesting issue in its own right in statistical 
mechanics, is of great interest in the holographic context. The typical states that dominate 
the ensemble have gross features of black hole geometries in the following sense: in quantum 
gravity it is these states that are expected to have a sensible semi-classical background. One 
issue however, is that our discussion will be confined to situations where we have states 
preserving two or more supercharges; this class of states (or their dual operators) in four 
dimensional superconformal theories do not include the supersymmetric AdSs black holes. 
Indeed, the entropy of the ensembles we consider scales as N a with a < 2; this is insufficient 
in the AdS / CFT context to produce sufficient back- reaction to obtain a semi-classical geom- 
etry with a regular event horizon. Nevertheless, BPS operators preserving higher amounts of 
supersymmetry provide a useful toy model to understand aspects of emergent gravity [37]. 

We begin with a brief introduction to the class of partition functions we focus on in §2 
and describe how these are related to vector partitions in §3. In §4 we discuss aspects of the 
partition function at finite N and show how one can recover explicit formulae exploiting the 



oo 





n=l 



3 



combinatorial properties of vector partitions. In §5 we analyze the density of states of these 
partition functions, generalizing the results of Meinardus. We then address issues regarding 
the typical states in ensembles under consideration in §6 using the notion of limit curves for 
the partition sums. Finally, in §7 we discuss some of the implications our results have for 
the understanding of supersymmetric operators in field theories in the holographic context. 

2 BPS cohomology partition functions 

The class of field theories we are interested in are superconformal field theories with a 
minimum of four supercharges, 7 on IR D for 3 < D < 6. As usual, operators in these theories 
are related to states onRx S 15 " 1 . These theories have some supersymmetry generators Q 
and superconformal generators S. Standard radial quantization results in the hermiticity 
conditions — S (see [38] for a nice discussion of radial quantization in SCFTs). The 
Hamiltonian H (the dilatation operator), the R-symmetry generators Ri (for a generic M = 1 
SCFT there would be just one) and the spins Ji (which come from the rotational symmetry 
group of M D ), all label the operators of the theory and satisfy (schematically) 



and therefore in a unitary representation one has a BPS bound H > ^\ -Rj+X^ Ji- Operators 
which saturate this are referred to as BPS. We will call any remaining global symmetries of 
the theory Gi (these are not present in M = 4 SYM, but generic M = 1 SCFTs may have 
them). A Q-cohomology 8 consists of operators that are annihilated by some subset of the 
Q that are not Q exact. This is equivalent to being annihilated by Q and (by analogy 
with Hodge theory). Therefore from (2.1) the Q-cohomology and BPS operators are in one 
to one correspondence, and this will be referred to as the BPS cohomology. Note that the 
BPS cohomology includes both primaries (annihilated by all S) and descendants with the 
same i?-charges. 

We will be interested in the partition functions over the BPS cohomology of at least 
two supercharges at weak non-zero coupling, with chemical potentials for each of the R- 
symmetries and global symmetries. The Hamiltonian is determined in terms of the charges 
as given in (2.1); denoting the Hilbert space of these BPS operators as 1-Lbps-, the partition 
functions generically take the form 



We count here just the Poincare supersymmetries; these theories contain superconformal symmetries as 
well. 

8 If one is considering a cohomology with respect to more than one Q then it is necessary that the Qs 
commute or anti-commute amongst each other, to ensure the cohomologies are compatible. 




(2.1) 



Z = TrH BPS exp(-rj Ri - g-id- Q Ji) 



(2.2) 
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where Ti,gi,Q are some set of chemical potentials conjugate to the symmetry in question. 
For example, such partition sums have been computed for BPS operators preserving various 
fractions of supersymmetry in Af = 4 SYM in [13], and for |— BPS states in the M2-brane 
CFT and |-BPS in the M5-brane CFT in [31], as well as a number of Af = 1 SCFTs in 
D = 4 [17, 23] for mesonic operators and in [25] for baryonic operators. Note that all these 
calculations assume that counting operators in the classical BPS cohomology gives the same 
answer as counting in the quantum corrected cohomology. 

The D = 4 case is of particular interest. 9 The Af = 4 result, was actually computed 
using a trick. In this case §-BPS states are chiral operators and thus their cohomology is 
equivalent to the chiral ring. Properties of the chiral ring could then be used to deduce 
that the finite N partition function over such operators is governed by an effective multi- 
dimensional harmonic oscillator model. The concept of a chiral ring in D = 4 only requires 
the existence of Af = 1 supersymmetry. Recall that for Af = 1 field theories it can be 
argued [40] that the chiral ring is generated by the gauge invariant combinations of chiral 
superfields in the theory <3>j and the gluino superfield W a , modulo some constraints. The 
constraints are the F-term constraints d<p i W = 0, with W being the superpotential, 10 along 
with [$j, W a } = and {W a , Wp} = 0. We refer the reader to [40] for an excellent account of 
chiral rings in Af = 1 supersymmetric field theories. One might expect that the chiral ring 
can be exploited to write down partition functions for BPS operators in Af = 1 SCFTs (as 
was done in Af = 4); this would involve solving the F-term constraint which depending on 
the superpotential can be a complicated affair. The "plethystic program" [21] and gauged 
linear sigma model techniques [23] have been used to deduce partition functions for the chiral 
ring in a wide variety of Af = 1 SCFTs (toric quiver gauge theories), some examples of which 
we will discuss shortly. 

Rather than being concerned with the specific partition functions that occur in given 
theories, in this paper we will study certain partition functions which can be used as building 
blocks of the known examples discussed above. There are two basic kinds corresponding to 
BPS bosonic operators and BPS fermionic operators. 

Let us consider the partition function for BPS bosonic operators at finite N, which we 
will denote by Zk(fl] N) where /3 = (/3i, • • ■ (3k) is a set of k chemical potentials conjugate to 
Hi, the quantum numbers of the various conserved charges in the superconformal theory in 
question and N is the rank of the gauge group 11 . The archetypal example of such a partition 
function is 

OO 

3*(/3,j0= z r-s — T (2-3) 

- LJ - 1 — v exp — p ■ n 



ni,ri2,"-nfe>0 



9 Systematics of the representation theory of the M = 2 and J\f = 4 superconformal algebras were devel- 
oped in [39]. 

10 For TV = 4 SYM the constraint reads = 0. 

n For simplicity we will focus on SCFTs whose gauge group is U(N). 
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where the infinite product converges if \p\ < 1 and Re (/%) > 0, which is actually a generating 
function for Zk(f3; N): 

oo 

E k (f3,p) = J2M(3;N)p N . (2.4) 

N=0 

Observe that this is the grand canonical partition function for bosons in a /c-dimensional 
harmonic oscillator potential, where p is the fugacity (the chemical potential that keeps 
track of particle number N). The partition functions given above correspond for instance 
to |— BPS or |— BPS states in Af = 4 SYM when k = 1, 2 respectively, or to |— BPS states 
states in the M2-brane world- volume theory for k = 4 and to —BPS states in the M5-brane 
world volume theory for k = 2 (the (2, 0) SOFT in six dimensions). 

Certain BPS cohomologies can include fermionic operators too. The partition function in 
this case typically consist of factors coming from bosonic operators (via products of E k (f3,p)) 
as well as factors in the numerator of the form 12 

oo oo 

S{(AP)= II (l+Pexp(-/3-n)) = ^Z/(/3;iV)^. (2.5) 

ni,n2,.-- n fc>0 N=0 

An important example of a BPS cohomology which includes fermions is the |-BPS cohomol- 
ogy for J\f = 4 SYM. It is given by [13] 

00 FT fl +T,pC* T 2i+l 2<J+1 2r+l\ 



l,q,r- 



; (l-px»2#xf) (l-pxf +2 X^ +2 X 2r+2 ) 



= H 3 (/3,p) S 3 (/3,px 2 x^) J] S|(/3,pe^x 1 x 2 x 3 ) (2.6) 

s=±l 

where in the second line we have shown how to rewrite it in terms of our basic partition 
functions (2.3) and (2.5) defining x\ = e~ ft . 

J\f = 1 superconformal quiver field theories dual to AdSs xl 5 with X 5 being a Sasaki- 
Einstein manifold have more complicated partition functions because of the presence of 
mesonic and baryonic operators (see [25] for a recent discussion of the moduli space and [36] 
for an analysis of typical states in such field theories). However, it has been argued that 
such partition functions can be written as sums over baryon number, with the summands 
each looking like "weighted" versions of the partition function (2.3) (i.e., each factor in the 
infinite product is raised to some n-dependent factor, which also depends on the baryon 
number). Hence we will also consider the more general class of partition functions 

3?0M= II 71 ) a nv5 (2J) 

n u n 2 ,n k > (1 "Pexp(-/3 • n))«- 

where w n is some sequence labelled by n (such that the infinite product converges for \p\ < 1 
and Re(/3j) > 0). For example, the partition function for mesonic operators (zero baryon 



J Generically, we will denote the fermionic analogue to a bosonic quantity by adding a superscript /. 
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number) in the quiver field theories dual to AdSs x Y p ' q is given by (2.7) with k = 4 
and w n = S n .Q t o where Q = (p + q,p — q,—p,—p) is a charge vector characterising the 
Sasaki-Einstein manifold Y p ' q [17, 23] (one can generalize to L abc - in this case the charge 
vector is Q = (a, b, c, — (a + b + c))). Of course, the weighted partition functions can be 
used to build the other partition functions so far discussed; for example if w n = — 1 then 
H{0M = S£G8,-p). 

We will be dealing with partition functions whose large N limit exists for all values of 
the chemical potential (i.e., it is O (l)). 13 Indeed the strict large N limit of Z k (f3; N) exists 
and we will call it 3fc(/3). The large N limit can be obtained via the following general trick: 

3 fc (/3)=lim(l-p)H fc 0M (2.8) 

p-*i 

which can be proved as follows. Consider 

oo oo 1 

J2 a NP N = (Zk(P,N) -3*09)) P N = S fc 09,p) - —3*09). (2.9) 

N=0 N=0 P 

From the fact that a N — > as iV — ► oo, it is easy to show that 14 (1 — p) ^jv=o a NP N - > as 
p — > 1 and hence the result is established. In the case at hand 

oo , 

3*08) = , (2.10) 

- LJ - , 1 — exp — 3 ■ n 

If there are also fermionic factors in the partition function, one simply sets p = 1 in these, 
and thus we define 3{(/3) = S{(/3, l). 15 

To summarise, there are three basic objects we will focus on; the grand canonical partition 
sum Hfc(/3,p), the finite N partition sum 2k((3',N) and finally the large N partition sum 
3fe(/3). As explained above we are also interested in the fermionic generating functions 
S{(/3,p) and 2>[(f3] N), however since their analysis is so similar to the bosonic case we will 
only indicate the necessary differences at the end of each section as appropriate. Similar 
comments apply to the more general weighted partition functions. 

One partition function which motivated our analysis is that which encodes the counting 
of the chiral ring |— BPS operators in Af = 4 SYM discussed above (2.6). In particular, one 



13 Certain partition function exhibit large N phase transitions for small chemical potentials. We will not 
consider such partition functions. We also only consider the large N limit where we keep the chemical 
potentials pi fixed. As discussed in [13], scaling the chemical potential with N in the large N limit results 

in an interesting Bose-Einstein condensation. 

14 This relies on Abel's theorem on continuity up to the circle of convergence, see e.g. [41]. 

15 Note that lim N z[(/3; N) = using (2.8); however such a limit never occurs as fermionic partition 
functions are always accompanied by a bosonic one which takes care of the factor 1 — p leading to a finite 
answer for the large N limit. 
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would like to be able to use the spectral data contained in the partition sum to construct 
an effective model encapsulating the dynamics of these BPS operators. Of greater interest 
would be to analyze the partition function of ^— BPS operators (these are outside the 
scope of our analysis as such operators are only preserved by only one supercharge in our 
language). In this case one knows the free field theory answer [13], as well as the count 
of planar operators i.e., graviton states with dimensions of 0(1), at strong coupling [13] 
(using supergraviton representations) and at non-zero weak coupling [42] (using the one- 
loop dilatation operator). The analysis we undertake is geared towards applications in the 
holographic context; in particular, to address questions regarding which operators (or states) 
are expected to be dual to semi-classical geometries and perhaps to aid in the construction 
of micro-state geometries from the supergravity side along of the lines of [43] for M = 4 
SYM. We will revisit these issues in the discussion after a detailed analysis of the partition 
sums in the next few sections. 



3 Vector partitions 

The finite N partition functions Zk(fl] N) and the large N partition functions 3fc(/3), which 
we have introduced in the previous section, have a combinatorial interpretation. This is well 
known to mathematicians, see e.g., [35] which we will base our discussion on. Let P(n) 
denote the number of partitions of n = (ni,7i2, • • -nk), an ordered /c-tuple of non-negative 
integers not all zero. These are often referred to simply as vector partitions. As for integer 
partitions one does not account for the order of the parts of the partition; to take care of this 
one introduces a concept of an ordering. More precisely one counts the number of distinct 
ways to write n = j 1 + j 2 + ■ ■ • + j s subject to the ordering j r > j r+1 , where j r > j r+1 if and 
only if jl > H +1 where i is the least integer such that jj ^ j? + . Note that this ordering 
allows us to write the condition ri\, n 2 , ■ ■ -n^ > and n^O simply as n > 0. If the number 
of parts in the vector partition is restricted to be at most N then we denote the number of 
such partitions by P(n; N). It is then a basic fact that: 

MP) = $>(n)e-^ n (3.1) 

n>0 

Z h (P;N) = ^P(n,iV)e"^ n . (3.2) 

n>0 

For completeness we will give the argument for (3.2) (the arguments for the other cases 
proceed similarly). Consider the frequency representation of a partition of n; that is n = 
S;>o a j J wnere a, denotes the number of times j appears in the partition. Observe that 



S 



J2j>o a i gi yes the number of parts in the partition. It is then clear that: 

oo oo 

J2J2 P ^ N ) e ~ f3n P N = E E expC-^^-jaj)^^^ 

7V=0n>0 a=0{aj|j>0} j>0 

1 °° 1 1 

= r^nE ( pe j = -^—p n 1 _ pe-p-i 

y j>0 Oj=0 y j>0 ^ 

which establishes the result. Note that in the first line each sum over aj is from to 00 and 
a = N — J2j>o a i ( reca h we are looking at partitions whose number cannot exceed N), and 
in the second line the various geometric sums have been performed. 

We will be interested in the asymptotic density of states for these partition functions, 
and from the above this can be seen to be equivalent to the asymptotic number of vector 
partitions. 

The fermionic partition functions also have a combinatorical interpretation in terms of 
vector partitions. Let Q(n) be the number of partitions of n (as defined above) into distinct 
parts, where (0,0,...,0) may be a part. Similarly denote the number of vector partitions 
of n into N distinct parts (including 0) by Q(n; TV). It then turns out that the generating 
functions for these are given by our fermionic generating functions 



3{09) = ^g(n)e^ n (3.3) 

n>0 

Z[(/3;N) = ^Q(n,iV)e^ n (3.4) 



n>0 



4 Finite N generating functions 



We will now explore the finite N partition sums in some detail. Our aim is to utilize some of 
the combinatorial properties of these partitions to give an algorithmic method to compute 
these from the grand-canonical partition function. The motivation behind undertaking this 
exercise is to learn about the operators that are present in the finite iV theory. Of course, 
this is relevant in the context of holography only when we are interested in quantum effects, 16 
but the analysis reveals interesting results which should enable one to better understand the 
set of BPS operators as we will discuss. 

We first note that the k = 1 case is well known: 

N 1 

n=l 

16 Recall that in the AdS/CFT correspondence, string effects are suppressed by g s ~ 
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which can be deduced by a variety of methods (for instance, directly using the well known 
free fermion description of the system or using the q-binomial theorem [24]). 

We may now exploit a useful trick in extracting the finite N generating function from 
Sk(f3,p). Note that the logarithm of the grand canonical partition sum admits a simple 
Taylor series in p 

_ _ „m„— m/3-n / k \ 



m — ' m 

m>l,n>0 m=l \i=l 



A general result is that for any function of G(p) such that 



m=l 

then 

F t) 

1 m y 



z — ' m! 

m=0 

where F m = Y m (gi,g 2 , ■ ■ ■ g m ) are the Bell polynomials [35]. They 17 can be defined by 



d m e y d l y 

Y m (yi,y2, ■■ -Vm) = e v -^r where yi = ~dtf- 

The first few are: 

}o = l, Y i(9i) = 9i, Y 2 (gi,g 2 ) = g\ + g 2 , Y 3 (g u g 2 , g 3 ) = g\ + 3gig 2 + g 3 (4.6) 
It can be shown that they satisfy the following recurrence relation: 



F^+i(<7i, g 2 , •■ ■ ,9n+i) — n C m Y n _ m (gi, ■ ■ ■ g n _ m ) g m+ i. (4.7) 

m=Q 

Thus, in the case at hand G{p) = Ef.(/3,p) and hence, from equation (4.2) we can read off 

u/9) = (»»-i)! n (4 - 8) 



i=l 



which gives the following expression for the finite N partition function 

Z k ((3;N) = l^Y N ( gi ((3),g 2 ((3),--- ,g N ((3)) (4.9) 

which is a result known in the mathematical literature [35]; this provides a generalisation 
of (4.1) for k > 1. For low values of N this is easily computed using (4.5) or the recursion 
relation (4.7). 



17 Not surprisingly, the Bell polynomials are closely related to the notion of the plethystic exponential 
introduced in [17, 21]. 
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Similarly, one can work out the finite N partition function for fermionic and weighted 

vector partitions. The answers in both cases are given by (4.9) with g m ((3) replaced by one 
of the following: for the fermionic partition function Z[(f3; N) 

gU/3) = (-l) m+1 g m (f3) (4.10) 

whereas for the weighted partition function 2^(/3; N) 

gZ((3) = (m-l)\J2wne- m(3n . (4.11) 

n>0 



4.1 Factorization property of finite N generating functions 

The Bell polynomials are very useful to infer certain properties of the finite N partition 
functions directly. They allow one to prove the following 

Factorisation property: The partition function Zk((3; N) can be written as 

Zk{P]N) = (jl^ifcN)] V k {/3;N) (4.12) 

where Vk{/3] N) is a symmetric polynomial in Xi = e _/3i of total degree at most \ kN (N — 1). 
In particular, when k is even 

/ k \ 

V k ((3;N)= I J[xA + 1 (4.13) 

where the • • • stand for the non-constant terms of the polynomial of lower order. 18 
Proof: We first note that (4.1) implies 

k N 

nw=iiS- (414) 



Then the recursion relation (4.7) allows us to deduce one for Vk(/3; N): 

N 



1 N 

V k (/3; N +1) = y~] a m , N g±g 2 --- g m gN-m+i V k ((3; m) (4.15) 

fU fir, . . . n»r , 1 ' 



m=0 



where a mt N = ( N+1 ^ N _ m -j\ Yl n =m n '- ^ ne can use ^ ms ^° P rove that Vk((3; N) is a polynomial 
by induction. Assume that Vk{f3; m) is a polynomial for m < N and observe the base case 



18 In fact this result is known in the mathematics literature [44], a fact we were not initially aware of; we 
therefore provide our own proof of this result. Further, [44] also proves that the coefficients of the polynomial 
are non-negative. 
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Vk{f3; 1) = 1. It is clear that all the terms in (4.15) for which N — m+ 1 > m are polynomial 
in Xi since one can cancel each g n in the numerator with one in the denominator, thus leaving 
only a product of some of the g n in the denominator (recall l/g n are polynomials in Xj). At 
first glance the terms for which N — m+ 1 < m are not obviously polynomial, as one will have 
a factor of g 2 N _ m+1 in the numerator but only a gN-m+i in the denominator. For these terms 
argue as follows. Cancelling g\ g 2 ■ ■ ■ g m in the numerator with the corresponding factors in 
the denominator will leave 

(416) 

9m+l 9m+2 ' " " 9N+1 

Note that m + 1, m + 2, • • • , N + 1 is a sequence on iV — m + 1 consecutive integers and thus 
one of these must be divisible by TV — m + 1. Lets call this particular integer M so that 
M = q(N — m + 1) for some integer g. Then since 



g N -m+i (N-m)\ 



9m 



(M 



1 ! M l-x 

' 1=1 



iV-m+l 



.r 



+ • • ■ + [X. 



i=l 



(4.17) 

is a polynomial we have proved that all the terms in (4.15) with N — m + 1 < m are also 
polynomial. Therefore we deduce that Vk(f3; N + 1) is polynomial completing the proof by 
induction. 

One can use the above argument to extract the order of the polynomial Vk{f3, N). This 
depends on whether k is even or odd. In particular, note that: 



9N-m+l 



9m+l9m+2 ■ ■ ■ 9N+1 



(N — m)\ 



JV 



111 



|(iV+m)(iV-m+l) 



X j 



(4.18) 



(which holds for all < m < N). First consider when k is even. We will prove (4.13) 
by induction. The N = 1 case is trivial. Now assume, for induction, that (4.13) holds for 
Vk(f3,m) for all m < N. Using (4.18), (4.15) and the induction step gives 



Vk((3,N 



N + 



1 N 

— Y 

+ i ^ 



m=0 

N 



i(A r +m)(A r -m+l)+im(m-l) 



— E 

+ 1 ^ 

m=0 



n 

k 

n 



X j 



+ ••• + 1 



ijV(JV+l) 



+ 1 



ijV(iV+l) 



n 



•1 i 



+ ••• + 1 



(4.19) 



and thus we learn that (4.13) also holds for Vk(f3,N + 1). Hence by induction this proves 
(4.13) for all N. 
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The situation with k odd is tricker, since the leading order term in (4.13) tends to get 
cancelled. While the precise order seems to depend on particular values of k and N it remains 
true that Vk((3] N) grows at most as fast as (xj) 5 7V ( 7V_1 ) j n each of the x\. 

Remarks: The above result has an interesting implication. The partition function Zj,((5; N) 
is that of iV-bosons in a /c-dimensional harmonic oscillator. The factorisation theorem we 
have proved shows it is equal to the product of k distinct partition functions of iV-bosons 
in a 1-dimensional oscillator Z\((3i;N), times a partition function whose Hamiltonian has 
a spectrum which is bounded above and below which presumably provides the interaction 
term between the k systems. We discuss some implications of this result in §7. 



Fermions: The partition function Z[(f3; N) can be written as 

ZlifrN) = (l[Z 1 ^ i ;N) ) V f k ((3;N) (4.20) 



where V[((3; N) is a symmetric polynomial in Xi = e of total degree at most | k N (N — 1) 
with no constant term (for N > 1). In particular, when k is odd 

V[(f3,l) = 1, (4.21) 
V f k (/3;N) = fjjxij + --- + 0, N>2 (4.22) 

where the • • • stand for the non-constant terms of the polynomial of lower order (so this 
polynomial does not have a constant term for N > 2). 

To prove this one can use the same technique as we employed in the bosonic case, which 
we will only sketch. We find that V[((3; N) satisfies the recursion relation (4.15) with 
replaced by a^ N m = (— l) 7V ~ m a7v, m , and the proof that V[((3; N) is a polynomial goes through 
in the same way. To extract the order of the polynomial one uses (4.18), which for odd k 
will contribute a factor of (— l^ N ~ m to the leading order terms and will cancel with that in 
a Nm , thus giving the same leading order answer as the bosonic case. However, the constant 
terms will cancel due to the alternating sign in a Nm - 

Observe that the k — 1 case is simply 

N 1 

Z{(/3;N)=e-^ N ~^l[—^ (4.23) 

n=l 

a result which is easily derived by other means (e.g., one can use a method used to derive the 
g-binomial theorem). This result makes sense as ^N(N — 1) is the ground state energy for 
iV fermions in a harmonic oscillator potential. The fact that the constant term in V[((3; N) 
(for N > 1) is absent makes sense as this system has fermionic statistics. 
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5 Asymptotic density of states 



Consider the partition functions Zk((3; N). From a statistical mechanics point of view these 
encode the degeneracy of states carrying charges n, which we will denote by ^^(n). There- 
fore from the combinatorial interpretation of this partition function (3.2) we immediately 
see that the density of states is given by c4,7v( n ) = -P( n ? N). Thus 



Z fc (/3;iV) = ^4,7v(n)e- n/3 (5.1 

n>0 

which can be inverted to give 

4,iv(n) = ] ' / ^Z k ((3;N)e^ 



(5.2) 



for some arbitrary bi G R+. We will also define dfc(n) = liniAr^oo c4,7v(n). We will be 
interested in the asymptotic behaviour of dk,N( n ) an d c4(n) as rii — > oo. Of course, for 
k > 1 there is more than one way of doing this. We will focus on the case where we send all 
the Hi — > oo at the same rate. Other cases, such as keeping some subset of the n; fixed, are 
more complicated, and we will outline how these can be dealt with at the end. 

The integral for ^^(n) at large n can be calculated using the method of steepest descent, 
which requires choosing the contour (and thus bi) such that we pass through the dominant 19 
saddle point in the direction of steepest descent. Note that choosing the dominant saddle 
means that contributions from any other saddle points will be exponentially surpressed. If 
we let bi — * + we expect this integral to be dominated by the — * + + iO region of 
integration, since the singularities of Zj~(f3;N) come from those of g m {f3) defined in (4.8); 
the strongest of these occurs at = 0. Since we are sending the — > oo at the same rate, 
we expect that will will need to send the j3i — > + at the same rate. Indeed from (4.9) one 
can easily deduce that in this case the leading order behaviour is 

ZJ3-N) \ . (5.3) 

The saddle point is located at an extremum of the exponent of the integrand: 

n l + ^logZ k ((3;N) = (5.4) 

and thus, using (5.3), one gets $ ~ N/rii. This shows that sending the $ — > + at the same 
rate is consistent with sending the rij — > oo at the same rate. 

We will now present a systematic analysis of the asymptotic behaviour of the density of 
states by generalizing the classic results of Hardy- Ramanuj an and Meinardus [34]. 20 As a 



19 We will only encounter cases where there is a unique dominant saddle point. 
20 See [45] and [46] for other generalisations. 
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warm up we start on the single chemical potential case obtained by setting (3i = (3 V i. We 
then discuss the general case of unequal chemical potentials. We will present the results for 
the large N partition functions 3&(/3) and indicate some generalisations at the end. 



5.1 Equal chemical potentials: Meinardus generalized 

When the chemical potentials coincide, the large N partition sums of interest are as in 
(1.4). These are similar to the weighted partitions considered by Meinardus [34, 45] with an 
important difference. To see this let us start with the Mellin representation of the partition 
function (1.4) 



oo oo oo 1 

iog3*o0) = - Yl a - lo s i 1 - e ~ nP ) = EE- a « e ~ mn 

n=l n=l m=l /r: ri\ 

i /"7+ioo 

= — / dsT(s)as + l)D(s)f3- s . 



' 7— 100 

The first line is standard and in writing the second line we have used the Mellin representation 
of the exponential 

1 /»7+ioo 

e~ x = — / dsx- s T(s) (5.6) 

' 7— ioo 



2ni 



and re-summed the series in m to obtain the Riemann zeta sum ((s + 1). Finally, the 
summation over n leads to the Dirichlet series, D(s), as defined in (1.5). This is the point of 
departure from Meinardus; as we will see the Dirichlet series of interest have k > 1 simple 
poles. 21 These don't seem to have been fully considered in the literature. 22 One reason 
to ignore the subtlety is that generically the dominant contribution to the density of states 
comes from the rightmost pole (a fact used in the recent analysis of [36]). We will now present 
the effects of having the k simple poles, which is quite simple to implement algorithmically. 

So far we have not specified the precise contour of integration in (5.5). For convergence 
of the integral we need the contour to lie to the right of any singularity arising from the 
Dirichlet series (or the zeta function). This requires that 7 > k since the defining expression 
for D(s) converges for Re (s) > k. 

To proceed it is useful to obtain an integral expression for the Dirichlet series 



D(s) = — dtt^git), g(t) = J2a n t n = - — -1 (5.7) 

r (*) Jo ^ C 1 -e T 



21 This occurs because the a n of interest (1.4) are polynomials in n of order k—1, and thus D(s) is a linear 

combination of k zeta functions (each of which have a simple pole) with shifted arguments 

22 See however [45] for discussions of spectra in toroidally compactificd theories, where zeta-sums with 

multiple poles arise. 
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which allows us to obtain an expression for D(s) valid in the whole complex plane 

D{a) = — \—± / dt{-ty~ l g{t) (5.8) 

where we have used a standard notation for the Hankel contour. 23 To obtain (5.8) we used a 
standard trick to the one employed in the analytic continuation of the Gamma-function, or 
indeed the Riemann zeta function [41]. The integral representation (5.8) allows us to deduce 
that D(s) is meromorphic with simple poles at s — 1, 2, ■ ■ ■ , k. Let the residues of D(s) at 
these poles be Aj where j — 1, 2, • ■ ■ , k. One can show that for the case at hand 

(— lV +fc k B n 

Aj = ReS s=j D( S ) = _ Yl ^n 1+ n 2 +...+n k ,k-j JJ -y (5.9) 

^ ni,...,n fc >0 i=l l ' 

where B n are the Bernoulli numbers, and thus note in particular that A/~ = l/(k — 1)!. 
However, we can proceed with the calculation for general Aj. 

We now return to (5.5), where we shift the contour at Re (s) = 7 to the left, across the 
poles of the integrand, to Re (s) = — e with < e < 1. Therefore log3fc receives contributions 
from the poles of the integrand in the region — e < Re (s) < 7 (which include those of D(s)) 
and we find 

3*09) = exp ^£ A ^j)aj + l) + D , (Q) _ m lQg/3 j {1 + a m y (5 1Q) 

as ->• + . 

We are now in a position to extract the density of states. First, observe that 3k(P) is 
actually a generating function for 



ni+ri2H hn k ,n d fc (n) (5.11) 



n>0 



as can be seem by setting — f3 in (5.1). Therefore we will derive an asymptotic formula 
for Cfc(n) as n — > 00. To derive this we use the integral expression (5.2) 



Ck(n) = — / d{33 k (/3)e n ? (5.12) 
valid for any b > 0. To obtain the saddle point evaluation of (5.12) consider 

k 

Gk((3) = J2^k+ n ( 3 - ( 5 - 13 ) 
.7=1 3 p 



23 The Hankel contour starts just above the real axis at +00, encircles the origin counter-clockwise, and 
ends up just below the real axis at +00. 
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Then, for (3 — > + (5.10) tells us that e Gk ^ provides a good approximation to the integrand 
in (5.12), with 

C j -> ; r(./ • 0a./ • 1). (5.14) 

We will choose the integration contour to pass through the saddle point of G k (P) with 
largest value of G k , along the direction of the steepest descent. Let this extremum be j3 n , so 
G' k ((3 n ) = 0, and thus we set b = (3 n . We have 

t5l = « (5-15) 

3=1 Pn 

and thus /3~ l is the largest positive root of this polynomial of order k + 1. This shows that 
n — > oo is equivalent to (3 n — > 0. More precisely 

/5„= (l + 0(n-^t)) (5.16) 

and it is useful to record 

cJ(W = £<^ = ^4^(i + e> (»-*)). (s.iT) 

Note that since G k (/3 n ) > the direction of steepest descent is in the imaginary direction. 
Thus, now we change variables to ir = \/G k ((3 n ) ((3 — (3 n ). This implies that as n — > oo (for 
fixed r) we have 24 

GM = G k (f3 n ) - ~r 2 + <D{p k J 2 ). (5.18) 

It is also useful to note that log (3 = log j3 n + O \J3nj. Putting all this together and we 
obtain 

i/3 -D(0) e G M+ D'(0) (l + (^)) 



c k ln 



/ , - f7WT7T -, exp G k ((3 n ) - D(0) \og(3 n - -r 2 + O (f3 k J 2 ) X (1 + O (/£)) 



2ttGZ(&) 

$^ eGkWnh-D'fp) (i + (n-^)) (5.19) 



(1-20(0)) 
C£ 2 < fc + 1) 2O(0)- fc -2 

n 



V^vr^ + l) 



where, in order to perform the integral, the limits of r have been replaced 25 by ±oo. The 
last equality follows from using (5.16) and (5.17); note that the pre-factor multiplying the 
exponential is the same as Meinardus' result [34] as it comes from the rightmost pole. This 



24 To derive this notice that G^\f3 n ) = O (/3~ p - k ). 

25 The error in doing this is exponentially suppressed since tt G'l(f3 n ) — > oo in the limit we are working. 
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is the desired asymptotic formula. Note, as we will show shortly, that we can obtain (3 n and 
Gk{j3 n ) as a function of n to the desired order using the defining polynomial for (3 n - however, 
only the positive powers of n in the expansion of Gk{P n ) are meaningful in our asymptotic 
formula (and there are only a finite number of these). Also note that to leading order 

logc fc (n)~G fc (/3 n )~ + n& [k({k + l)\^ (5.20) 

which agrees with [24]. 

Finally, observe that the above derivation is valid for any Dirichlet series D(s) with k 
simple poles which converges for Re (s) > k with residues - in fact the number of poles 
and the region of convergence need not be related. 26 Thus specialising our formula to the 
case of a single simple pole at s — k, one can check that it correctly reduces to Meinardus' 
theorem. 



Explicit formulae: We have shown how the asymptotics of Cfc(n) may be determined up 
to some known function Gk{f3 n ). We will now address the explicit calculation of Gk{f3 n )- 
First consider k = 2. In this case satisfies a cubic and thus we can get an analytic 
expression for it in terms of n. Explicitly we have 



n 







We want the largest positive root to the above equation. Solving the cubic gives 



(5.21) 



where 



V 



n 
2G> 



3Co 



Ci 



3Co 



3C 2 



(5.22) 



\ 



1/3 



n 
2C? 



3C, 



3C 9 



(5.23) 



/ 



From this it is a simple matter of computing the expansions of fl n , G((3 n ) for large n which 
gives the explicit formula 



1 D(0) 

3 



C2(n 



-713 



i(D(0)-2) 



exp 6 -Cl'\^ + 



c. 



1/3 



n 



1/3 



ci 

QC, 



(5.24) 



26 0ne can take a Dirichlet series which converges for Re (s) > k and has an arbitrary number of poles in 
the region < Re (s) < k; each of these poles then would contribute to the integrand in the saddle point 
evaluation. 
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Observe that this illustrates the fact that for the leading asymptotics of Ck(n) one needs more 
that just the contribution from the rightmost pole of D(s). To see this, apply our formula 
(5.24) to a partition function whose Dirichlet series has C\ = and call the density of states 
c 2 (n); it is then clear that lim^oo c 2 (n)/c 2 (n) 7^ 1, i.e., it is not true that c 2 (n) ~ c 2 (n). 

One can obtain explicit formulae for higher values of k by a similar method. This is 
despite the fact there is no general formula for the root of a polynomial of order higher than 
4. To get round this one can use the Lagrange inversion formula which tells one how to 
invert a Taylor series. Define A = and (3 n = y(X), so y — > as A — > 0. From (5.15) 



A 



y 



k-j 



fc+i 



(5.25) 



The Lagrange inversion formula then gives 



1=1 



1 

U 



ji-i 



(5.26) 



y=0 



Notice that if I is a multiple of k + 1, so / = (k+l)m for some integer m, (f)(y) 1 is a polynomial 
of order (k — l)m and thus h\ = (since I — 1 — (k + l)m > (k — l)m). Therefore 



l>0, 1^0 mod (fe+1) 



(5.27) 



from which it is a straightforward matter to extract the large n expansion of G{j3 n ) for any 
k. This shows how the computation of the positive powers of n in the large n expansion 
of G(j3 n ) is purely algorithmic. To illustrate this we give the k — 3,4 cases explicitly. For 
k = 3: 



0n 



4Cf 3 

-«4 



G 



2 _2 

T n 4 



A CI 

r^ 2 + 

2C| 



1 



Co 



24 C, 



ci 
8 a 



_ 3 

n 4 



C 1 



ci 



On 



1 

n 4 



On 



(5.21 



(5.29) 
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and for k = 4 we get: 

5C| 5Q V &C 4/ 



6 



^(j5 \ 5 C4 25 C 4 



5Q i C3 3 1 / C| \ 2 1 / C*2 C3 C| 



4 3C| 2C| V 5 W d V 5 °4 25 G 4 



5Ca V ' 2Ca 2 12C 2 



+ ^ - ^ " 7777^2 + K 5 • ( 5 - 31 ) 



Fermions: The partition function "5{(P) can be represented by an integral of the form (5.5) 
with C(s + 1) — > Ca(s + 1) where Ca(s) = 2~^^=i( — t) n+1 /n s is the alternating zeta function. 
Since (a(s + 1) = (1 — 2~ s )((s + 1) we see that (a has no poles. Therefore the integrand of 
(5.5) will have a simple pole at s = (as opposed to a double one as in the bosonic case). 
This leads to the (3 — > + asymptotics being slightly modified: 



(5.32) 



and thus 



4(n) = \l 27rGk „ m e G ^ +D ^2 {1 + m) (5.33) 



where Gk(/3) is given by (5.13) with the coefficients Cj = AjT(j + 1) Ca(j + !)• Again 
Gk{f3 n ) = corresponds to the dominant saddle point. Observe that the leading order 
behaviour 

\ogc{(n) ~ G h {/3 n ) ~ fl + ^ \ k U(k + l)]^r (5.34) 

is the same (up to numerical factors) as in the bosonic case. Also note that one can use the 
explicit formulae developed above with the appropriate choice of Cj. 



5.2 Non-equal chemical potentials: multivariable generating functions 

Now we will consider the asymptotics of the large N generating functions 3fc(/3), for non- 
equal chemical potentials fy, as well as their associated density of states dk(n). It turns out 
that these generating functions can be written as products of the following functions: 

00 1 

Fk{P) = -, ; a c- (5.35) 

AJ - 1 — exp(— 3 ■ n) 
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More precisely, 

k 

3*(/3)=n n F i(Pii>Pi*>- (5.36) 

j=l ii<i2<—<ij 

Thus it is convenient to concentrate on the Fk{(3) as they provide building blocks for the 
generating functions we are actually interested in. Convergence of the infinite product is 
guaranteed by taking Re (/%) > 0. We will now derive the behaviour of this function as 
j3i — ► + + iO. Of course this limit can be taken in a number of ways: we will be interested 
in when all the approach zero at the same rate, so /3j = 7^ with (3 — > + and 7$ fixed. 
Taking the logarithm and expanding as before leads to 

logF fc (/3)= ^ —— = — dsF(s)C(s+l)C(s;(3) (5.37) 

where the second equality follows from using the inverse Mellin transform of the Gamma 
function (5.6) to replace the exponentials, and we have defined the following generalisation 
of the Riemann zeta function: 

00 1 

£ (^)7 (5.38) 

which converges for Re (s) > (so -y > k). 

We will need the full analytic structure of ((s, (3) as a complex function of s. Fortunately 
this can be worked out in a manner analogous to the zeta function, or the Dirichlet series 
discussed in §5.1. This involves using the integral representation of the Gamma function to 
derive: 

k 

1 /*oo 2 I 

which again is only valid for Re (s) > k and then one can analytically continue this expression 
using the Hankel contour to get 

C(«, f3) = dt (-ty-i n ^T3i - - r ( x - s ^ ®- ( 5 - 4 °) 

Observe that I(s, (3) is an entire function of s. It vanishes for s — k + 1, A; + 2, • ■ • . Therefore 
the apparent singularities coming from T(l — s) at these points are in fact removable and 
therefore ((s,f3) is analytic at these points. This leaves s — 1, 2 ■ • • , k as potential poles, 
and indeed it is easy to see I(s,(3) will be non-vanishing at these points, and thus they are 
simple poles. Therefore the above expression provides the analytic continuation of £(s, (3) 
to the whole complex s plane, resulting in a meromorphic function with simple poles at 
s = 1, 2 • • • , k. The residues of these poles are 

Re Ss=3 C(a, (3) = JJ^JJ ^ ^ ( 5 - 41 ) 
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for 1 < j < k where 

k 

I( 3n ) = Res t=0 (-ty- 1 ] 



-IV'- 1 



i=l 



rii=l li ni 



s <* 

712, ■■• ,1fe>0 



+ri2H \-n k ,k-j 



TT— 7 



The second equality follows from using z/(e z 
Bernouilli numbers B n . Note in particular, 



■i) 



fc-i 



(5.42) 

J2^=o B n z n /n\ which defines the 



(5.43) 



From the defining expression for £(s, /3) it is easy to see that ((s, (3) = 0{(3 s ) provided 
we avoid the poles. This implies that for all < e < 1 



dsT(s)as + l)as,(3) = 0^ 



(5.44) 



as (3 — »• + , since the integrand has no poles in the —1 < Re (s) < region of the complex 
s-plane. Therefore the asymptotic behaviour of Fk((3) as (3 — »• + is simply given by the sum 
of the residues of the poles of the integrand in (5.37) in the region — e < Re (s) < k with the 
error in this estimate given by (5. 44). 27 The result of this analysis is: 



F k {(3) = f3 T ^ exp 



E CO' + 1) IU, 7) P~ j + C'(0, 7) 

.3=1 



;i + (5.45) 



as (3 — > + . Note that to derive this we have used the following logarithmic property of our 
zeta function (5.40) 

C / (0,/3)=C / (0, 7 )-/(0,7)log/3. (5.46) 

Using these results it is a straightforward matter to deduce the asymptotics of 3fe(/3) 
using (5.36). For example, for k = 1: 



3 1 (P)=F 1 (P) = \l— exp 




L 

2tt 



7[~ 



Qf3 



(5.47) 



which is the well known result of Hardy and Ramanujan. The k = 2 case is new and gives: 
3 2 (/3) = F 1 (f3 1 )F 2 (p 2 )F 2 ((3) (5.48) 



/5* V 3 72 371 



exp 



r c(3) 


/ 7i + 72\ 


.7172/3 2 


V 2 7i72 / 



C'(0,7) 



;i + o{(3*)) 



2tt 

where /3 = (ft, /3 2 ) = 

27 Note that the contribution from the paths joining the contours Re (s) = 7 and Re (s) = — e vanish due 
to the exponentially decaying behaviour of T(s) in the imaginary direction. 
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For higher k it is convenient to define the following zeta-function: 



k 

T ^) = Et^ = E E C{s, fa, (5.49) 

which also converges for Re (s) > k. We may exploit the analytic continuation of ((s, (3) to 
deduce that of Y(s, /3) in the complex s-plane, which is given by 

T(s,f3) = -T(l-s)l(s,(3) , (5.50) 

where 

k 

I(s,(3) = J2 I&PinPi*,--- ,fa)- (5-51) 

3=1 ii<i2<—<ij 

This shows that Y(s, (3) also has simple poles at s = 1, 2, • • • , /c. 28 We can therefore write 
down the asymptotic formula for 3fc(/3) in a fashion analogous to (5.45) for Fk(f3); one has 



3 fe (/3) = exp 



' k 

^(-iy +1 CO' + l)m/3) + T'(0,/3) 
.5=1 



<W) • (5-52) 



This asymptotic formula is one of the main results of this section. 

Now we are in a position to turn to the asymptotics of the density of states c4(n) of 
3fc(/3). Thus we wish to evaluate 

4(n) = ] / ^ 3*09) e n/3 (5.53) 

i=1 Jbi-m 1711 

in the limit rii — * oo (at the same rate), which we will do by the saddle point method. The 
asymptotic formula (5.52) tells us that e Gk ^\ where 

k 

G k ((3) = n-[3 + J2 (-1) J+1 C(j + 1) Ah 0) , (5-54) 

5=1 

provides a good approximation to the integrand for small We will choose the path of 
integration to go through a saddle point of Gk defined by diGk = such that Gk is largest. 
We thus choose the 6, appearing in (5.53) to be this particular saddle point, which defines 
the hi as functions of the As we will see shortly, the Hessian Hij = didjGk evaluated at 
the dominant saddle point is positive definite (at least in the limit of interest), and therefore 
the path of steepest descent is in the imaginary direction. Now, observe that Hij(b) is a 
positive definite symmetric real matrix and thus must have real positive eigenvalues A* such 



28 0f course, log3fc(/3) can be represented as integral (5.37) with ((s,/3) — > T(s,/3), and we could have 
proceeded from here from the beginning. For clarity, though, we chose to analyse the Fk{(3) individually. 
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that Hij(b)ej = X e\ where e*- is a set of orthonormal eigenvectors. This allows us to change 
variables from $ to Tj defined by $ — bi = iTjXj so that 

G k ((3) = G k (b)-^r l r l + --- (5.55) 

where • ■ • denotes terms in the Taylor expansion of higher order. We can bound these higher 
orders terms as follows. Note that from diGk{b) = 0, which is a polynomial of total degree 
k + 1 in the 1/bi, the limit rii — > oo (all at the same rate) occurs when bi — > (at the same 
rate). In fact for small $ 

G fc (/3)~n-/3 + ^!±^ (5.56) 

1 Lj=l Pj 

from which one can compute all derivatives d^d^ ■ ■ -di p G k (b) in the limit bi — > 0. Using 
diG k (b) = this leads to: 

6< , I bj — 5.57 

which shows explicitly that if we send all the rij — > oo at the same rate then 6j — > at the 
same rate. Observe that the Hessian in this limit is 

from which it follows that 

det[H l3 ((3)] ~ (k + 1) / (fc + 1 ^ +2 • (5.59) 



Therefore Hij(b) is positive definite (at least in this limit) as promised earlier. We let b = be 
where b — > and c is fixed as rii — >• oo; in particular, 6 = O (\b\). Then it is easy to see that: 

d tl d l2 -A p G k (b) = 0{\b\- h - p ). (5.60) 

Thus from Hy(b) = O (|b| _A;_2 ) we see that the eigenvalues A J = O (|b|~ fc ~ 2 ) and hence for 
fixed Ti we have # - 6< = O (|bp+ 2 )/ 2 ) , so 

ftA ■ • (Ai - 6 n) • • • (A P - U = (|&| (p - 2)fc/2 ) (5.61) 

and therefore for p > 2 these terms are of order (9 (|6| fe / 2 ); this gives us a bound on the 
terms denoted by the ellipses in (5.55). One can also argue that for fixed Ti 

T'(0,/3) = T'(0,6) + 0(|b| fc/2 ) (5.62) 



21 



using (5.46). Therefore we have established that for fixed Tj 

log3*(/3) + n • (3 = G k (b) - ^r in + T'(0, b) + O (|b| £ ) (5.63) 

in the — » oo limit which is of interest to us. Taking the Jacobian for the change of variables 
into account allows us to deduce 

4(n) = \ =exp(G fc (b) + T'(0,6)) (1 + 0{\b\ k ' 2 )) (5.64) 

v /{27r) k det[H ij (b)\ 

C(fc+1)fc+1 m ) exp(G fe (6) + T'(0,b)) (5.65) 

where the second line follows from using (5.59) and (5.57). This asymptotic formula is the 
main result of this section. Note that to leading order 

k 

fogd k (n)~G h (b)~(k + l)t(k + l)T& flnf 1 (5.66) 

i=i 

which agrees with the result obtained in [24]. Also note that in the k = 1 case (5.84) 
simplifies to the classic result of Hardy and Ramanujan (1.3). 



5.2.1 Generalizations 

Finite iV: Now a word on the finite N partition functions Zk((3, N). The k = I case was 
derived in [47] (see [21] for application to superconformal theories), and our discussion will 
include this as a special case. There are a number of possible methods one can contemplate 
to extract the density of states dfc,jv(n). Logically, the most straightforward would be to 
determine the small $ asymptotics of Zk((3,N), and then use the saddle point method to 
deduce the asymptotic density of states, just as we did above in the large iV case. However, 
the method we used above to calculate the small $ asymptotics does not obviously lend 
itself to Zk(/3,N) as we do not have an (infinite) product representation of it, the best we 
could do was write it in terms of Bell polynomials (4.9). A more indirect method, would 
be to calculate the small $ asymptotics of Ek((3,p), for fixed p, and use this to deduce the 
small Pi asymptotics of Z({3, N) using 

= M£ (5 - 67) 

The advantage of this, is that the method applied to the large case 3fe(/3) generalises 
straightforwardly to Hfc(/3,p). In fact, following what we did in the large case, note that 

k 

2 fe (/3,p) = r ^n 11 F i(Ai,A 2 ,--" ,Pij,p) (5.68) 

■P j=l ii<i 2 <---<ij 
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where 

which can be written as: 
log F k ((3,p) 



F h {P,p) 



E 

m,ni>l 



n 



ni,ri2---n fe >l 



1 

1 — p exp(— (3 ■ n) 



-rafiirii 



m 



1 

2m 



7+ioo 



dsT(s)U s+1 (p)as-[3). 



(5.69) 



(5.70) 



7— JOO 



Note that this is of the same form as for log F k (P) (5.37), with ((s + 1) — ► Li s+ i(p). Recall 
that ((s + 1) has a simple pole at s — 0. In contrast, the polylogarithm function Li s+ i(p) 
is an entire function of s for |p| < 1 (since the defining sum converges absolutely for these 
ranges). Therefore the pole structure of the integrand in (5.70) is simpler than in (5.37); 
the only change being one has a simple rather than double pole at s = 0. One can go on to 
derive an analogous asymptotic formula for F k (P,p), which gives 

k 



F k ((3,p) = exp 



U j+1 (p)lV,<y)j3- 



j=Q 



(5.71) 



as P — > + . This can then be used to deduce the asymptotics of Ek(P,p) using (5.68). 
Now we turn to the evaluation of c4,7v( n ) given by (5.2). Using (5.67) and assuming we 
can swap the order of the Pi integrations and the ^-differentiations allows one to deduce the 
asymptotics of ^^(n) for large n via saddle point integration in /?, as before resulting in: 

/ / _h_ X 1 / 2 \ 

Li fc+ i(p)fe+i 



1 d 



N 



4,Ar(n) ~ — 



N\ dp 



N 



V 



(2n) k (k + l)(l-p) 2 U 



k + 2 
k + 1 



exp[G k (b,p)] 



J 



p=0 



where 



G k (P,p) = n-P 



k 

E 

3=0 



iy +1 u j+1 ( P )i(j,p). 



(5.72) 



(5.73) 



As before b is given by the location of the dominant saddle point diG k {P,p) = and thus 
depends on p as well as n; to leading order b is given by (5.57) with ((k + 1) —> Li k +i(p)- It 
would be interesting to determine whether the steps leading to (5.72) are indeed valid and 
if so determine more explicitly its N dependence. 



Fermions: Let us mention what happens in the case of a fermionic partition function. 
Thus define 

H(P)= II l + exp(-/3-n) (5.74) 

ni,ri2,...n k >l 

which can then be shown to satisfy (5.37) with £(s + 1) — > Ca(s + 1) = (1 — 2~ s )£(s + 1). 
Therefore the pole at s = is simple in this case and one gets 

" k 

F/(/3)=exp J2(-l) j+1 Uj + l)nj,l)P- j + (5-75) 
.3=0 
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as f3 — > + . This allows one to derive 

/ -a. \ 1/2 

4(n) ~ Uk + l)^ exp[G f k (b)\ (5.76) 

V(2^(A; + l)n- =1 ^ + 7 

where 

k 

G{((3) = n-(3 + J2 UJ + l)I(j,/3) (5-77) 

j=0 

and the dominant saddle diG{{b) = defines b (as a function of n), whose leading asymp- 
totics are given by (5.57) with — > (^(fc+l). As in the bosonic case G{{b) = O (\b\~ k ) 
as 6 — > + , and hence 

logtZ{(n)~GfJ(b)~(fc + l)a(A; + l)^i ^nf 1 (5.78) 

i=i 

and thus to leading order behaves in a similar way to the bosonic partition functions. 

Weighted: In this case define 

FJf(m = TT i— — - (5.79) 

ni,ri2,--- ,rifc>l 

which can be shown to satisfy (5.37) with C( s ; /3) — ► C"'( s )/3) where 

and we will also consider Y w (s, /3) defined as in (5.49) with £(s, /3) — > £ w (s, The domain 
of convergence of these generalised Dirichlet series depends on u> n ; let us suppose ( w (s,(3) 
converges for Re (s) > a. One can represent this by an integral, valid for Re (s) > a as we 
did in the simpler cases 

1 r°° 

C{s ^ ) = T(s) dtirl 9°®> 9 w (t)= ^e-^K (5.81) 

\ S > ^° m,n 2 ,- ,n k >i 

The analytic continuation can then be performed using the Hankel contour: 

i Ho + ) 

r( s ,/3) = -r(i- s )/-( s ,/3), r>{s,p) = — / dt (-*)'- V(*)- (5.82) 

The contour integral I w (s, (3) defines an entire function of s. Thus the poles of ( w (s, (3) must 
come from the Gamma function, and hence can only occur at the positive integers. Exactly 
which integers will depend on the analytic structure of g w (t) as t = 0. If g w (t) has a pole 
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of order a, which is what we will assume, then ( w (s,f3) will have at most 29 simple poles at 
s — 1, 2, ■ ■ • , a since I w (s, (3) — for all integer s greater than a. Then one can repeat the 
calculations to get the asymptotics of F£(f3) as (3 — > + : 



F™({3) =/? /tu (°^exp 



.i=i 



(l + 0(fa)). (5.83) 



which then can be used to get the asymptotics of 3™ (/3) in the same way as (5.36). From 
this one can derive: 

1 



n 



(27r)*det[tfg(b)] 



exp (G£(6) + Tr"(0, 6)) + 



(5.84) 



as n -> cx), where G£(/3) are as in (5.54) with Z(j, /3) -> J w (j, /3), and ^G^(b) = defines b 
(as a function of n). Note that T w (s, (3) is the weighted generalization of (5.51). The leading 
asymptotic behaviour of logrf^(n) is given by G™(b), which in turn will be determined by 
the j = a term. However, given in general one does not know the exact form of I w (a,(3) 
(and thus b), we cannot be any more explicit at this stage without specifying w n . Of course, 
given a specific w n one can easily apply the above formalism to get explicit answers. 



Other asymptotics: We have only focused on asymptotics where the rti are all large and 
comparable. It is interesting to ask what occurs if one sends only some of the rii to infinity 
and keep the others fixed. It would then appear that one cannot use the saddle point method 
for the integrations in fa directions which correspond to the fixed n^. Thus, at best, one 
can only perform a saddle point integration in the remaining directions, which require the 
knowledge of the fa — > + asymptotics of the generating function in those directions alone 
(i.e., the other fa fixed). We will now outline how these may be worked out using the 
same techniques as before. The crucial step is to derive a different integral representation 
for log Fk(f3) which is adapted to the asymptotics we require. Thus suppose we want the 
asymptotics as fa — * + for i = 1, ■ • • , p while keeping fa for i — p + 1, ■ ■ • , k fixed. Expand 
logi 7 ) c (/3) as in (5.37) and use the inverse Mellin transform of the Gamma function to replace 
the exponentials e~ ni/3i7n only for i = 1, ■ • • ,p. One can then resum the remaining e~ nil3im 
for i — p + 1, • • • , k obtaining: 

logF fc (/3)= / T(s)L(s + l ] f3)C(s,(3') (5.85) 

where we denote (3 = (fa, ■ ■ ■ , (3 P ) and (3 = (fa+i, ■ ■ ■ , fa) and we have defined the function 

m=l j=p+l 

29 We say "at most" as it could occur that I w (s,(3) vanishes at some of these points thereby removing the 
singularity. 
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Note that, just as for Li s (p) for \p\ < 1, the series defining L(s;j3) converges absolutely 
for all complex s and thus defines an entire function in the complex s-plane. Therefore, as 
before, we may work out the asymptotics for small (3' resulting in: 



F k (/3) = exp 



i=i 



0{F)) 



(5.87) 



where (3' = f3~f' and (3 — > + . From this, using (5.36), one can deduce the asymptotics 
for 3fc(/3) itself in this particular limit. Now let us turn to the question of interest: the 



asymptotics of c4(n) in the limit of large n' 
Note that one can write: 



d fc (n) 



-A- l d n > 



n 



j=i 



[ni, 



n p ) with fi 



n 



■p+i) 



d(3i 



hi — Z7T 



27TZ 



: MP) e n, (3 ' 



rih) fixed. 



(5.88) 



where Xa 



If one assumes that the large n' limit and the Xj differentiations can 



be swapped, then one can perform a saddle point integration in the (3' direction (as this 
localises around (3' — 0); this then requires the asymptotics of 3fc(/3) for small j3' which we 
have worked out above. Carrying out such steps gives the required asymptotics of dk(n) in 
terms of Xj derivatives which we will not write down. It would be interesting to determine 
whether this procedure is in fact valid. 



6 Limit curves and typical operators 

Let us now consider the partition functions 2^(3, N) and 3k{f3) from a more statistical point 
of view. We would like to infer from these partition sums the class of "typical operators" 
which dominate the ensemble under consideration. In general from all the states in the 
Hilbert space H-bps we will find that one can identify a limit curve in the charge space on 
which the typical operators lie. Our analysis follows the statistical treatment of [48, 49], 
which has previously been used in the discussion of |— BPS states in M = 4 SYM in [50] 
and more recently for quiver M = 1 theories in [36]. 

To begin with let us introduce a probability distribution on the set of quantum numbers 
n by: 

so Y2n>oP k ( n i N) = 1- We may now define the expectation value and variance of n in a 
standard fashion 

(ni) = —log Z k {/3,N) (6.2) 
varH = (n?)-(ni) 2 =Qp]ogZ h {P,N). (6.3) 



29 



We will be now interested in features of this distribution in the limit of large (rij) all of 
the same order. Physically the $ have the interpretation of the inverse temperature of a 
canonical ensemble and thus we expect the only way to achieve (nA — > oo is to send /3j — > 
(at the same rate). However this temperature does not have any meaning intrinsic to the 
theory which the states belong to, it simply sets the mean "energy". In any case, as $ — > 
we can use our leading order result for Z k (f3, N) (5.3) to obtain 

/ v N r , N 

(n*) ~ -r, var[n,] ~ -j (6.4) 

Pi Pi 

and therefore 



var 



n, 



-=. (6.5) 

(ni) V^V 

Note that these expressions are valid for all N in the /3, — > limit. Therefore we see that 
in the large N limit these — > distributions are sharply peaked around the mean. This 
allows one to introduce the concept of a limit curve [48, 49] which we now discuss. Let us 
work with the strict large N limit partition function 3fc(/3) since as we have just argued this 
is the regime where the distributions localise around the mean the most. One finds the exact 
answers 

r i \- n| 2C(fc + l) 

var nJ = > „ . ~ ? . 6.7) 

^4sinh 2 (^) 

where the asymptotics for $ — >• can be deduce from our asymptotics for Fk(/3) (5.45), 
noting that log3fe(/3) ~ log-^fc(/3)- It thus follows that 

which indeed tends to zero as predicted above from the large N limit of the finite iV answer, 
although this also tells us how fast it tends to zero. As just discussed in the j3i — > limit these 
distributions localise around the mean and due to the combinatorial interpretation in terms 
of vector partitions we may read off the so called completion numbers of the "partition" 
of (rii). For a vector partition of say v, the completion numbers r v (n) are defined by 
v = Xln>o rv ( n ) n ' *- e -' they give the number of times n appears in the given partition of v. 
Thus, we may formally define the completion numbers of (nA, despite this not necessarily 
being an integer, by 

^)( n ) = ^r T - ( 6 - 9 ) 

Now, given a set of completion numbers one may define the following curve, or surface, in 
jgfc+i gj ven by , t 2 , • • • t k , 4>{t)) where 



5>v( n )- ( 6 - 10 ) 



n>t 
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Observe that f£°dti <f>i(ti) = Vi where we have defined <f>i(ti) = 0(t) ij=Oj yi. It is clear that 
4>(t) is not a smooth, or even continuous function for vector partitions of v. However, the 
curve one obtains using the "completion numbers" of (ni) (6.9), is a smooth function of t. 

In the limit of interest, (n*) — > oo with — ► 0, one may replace the sum defining 0(t) 

with an integral. To do this explicitly, since J °° dij</>j(ij) = we first need to perform a 
rescaling in order that the limit may be taken. Thus, define 



n 7=1 & , ( h tk 



so that 

(^)a nj=i 



(6.11) 



d*i = T7r~~TT - 1 (6-12) 
as Pi — > 0. Then, defining yi = PiUi, we see that: 

*> - ^gM - t n ^^ry - cornj^ (613) 

which defines the limit curve C(t). Hence for the bosonic partition sum (2.10) 
Observe that J °° (it, Cj(ij) = 1 which was guaranteed by the construction. 



Fermions: For the applications we have in mind it is also useful to obtain the fermionic 
counterpart of the limit curve C(t) defined above. From the fermionic partition function 

3{(/3)=n i+e " /3n ' ( 6 - is ) 



n>0 



one obtains 



vara 



/ \ \^ n i U{k + 1) . . 

<«<> = Eto-^j^. ("6) 

= V n l ~ 2C ^ + 1 > (617) 



n>0 



4 cosh 2 (^) PHlUfr 



where the asymptotics for Pi — »■ + are deduced from our asymptotics for F[{(3) (5.75), 
noting that log3{(/3) ~ log Fj[(0). Then 



k 



va(Hi). 



I/?] 72 (6.18) 
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which shows that in the $ — > limit this distribution is sharply peaked around the mean 
as in the bosonic case. Therefore in this limit the concept of a limit curve is sensible and 
repeating the steps above leads to the limit curve: 

Applications: Having obtained the behaviour of the limit curves for the basic partition 
functions of interest, we now turn to some applications of this concept. Observe that the 
bosonic k — 1 case corresponds to a curve in M 2 ; one way to interpret this is to think of it as 
describing the behaviour of Young Tableaux's in the limit described above (recall that the 
standard integer partitions are naturally associated with Young Tableaux). In particular, the 
limit curve defined in (6.14) for k — 1 is basically the typical Young tableaux in the ensemble. 
This is related to the concept of "typical" |-BPS operators in M = 4 SYM [50]. 30 From the 
dual supergravity in AdS 5 xS 5 description |-BPS states correspond to smooth geometries 
which are also specified by a curve in M 2 (the LLM plane). Indeed, using the limit curve 
as a boundary condition for the supergravity solutions leads to the so-called "hyperstar" (a 
singular geometry) intepreted as the effective geometry dual to generic |-BPS operators of 
fixed charge [50]. 

An interesting question is whether such a mapping occurs in the |— BPS and §-BPS 
sectors. The k = 2 bosonic case above gives the limit curve for |— BPS operators and 
consists of a 2-surface in M. 3 . The supergravity point of view is far less developed than in 
the |-BPS case, essentially due to the non-linearity of the resulting equations - see [52] for 
the supergravity analysis. However, it has still been argued that smooth geometries should 
exist, specified by a smooth 3-surface in a 4d Kahler space [53] (see also the recent analysis 
of [54]). This therefore raises a puzzle as the boundary conditions for these supergravity 
solutions consist of a surface of different dimensionality to the limit curve. 

The |-BPS case as we have been earlier includes contributions from fermionic partitions 
functions. In order to construct the limit curve for this observe the following fact. Suppose 
we have two partition functions 3i arid 32 which generate some class of vector partitions 
with associated limit curves C\ and CV The limit curve associated to 3i 32 is C\ + CV 
Armed with this and the above results we can try and compute the limit curve of the large 
N |-BPS partition function (2.6). The special case = I/X1X2X3 is straightforward as it 
can be written in terms of our basic partition functions: 



1 U<j Mpu Pj) Yli<j % (A, Pj) 

Thus the limit curve in this case is: 

C , ( 1 / 8 )(t) = 2C3(t) + 2C 3 / (t)-^ [<7 2 (M;) + C£(Mi)] • (6-21) 



*<3 
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This limit curve is also relevant for thermal ensembles in a free boson theory in two dimensions [51]. 
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It would be interesting to work out the limit curve for general chemical potentials, but we 
shall not pursue this here. The |-BPS limit curve above consists of a 3-surface in a M 4 and 
presumably including the extra chemical potential would only increase its dimensionality 
by one. The supergravity situation in this sector is even less understood [55], but it seems 
smooth solutions should be specified by 5d boundary conditions in some 6d Kahler space 
[53] and thus again it is not clear how this relates to the limit curve we have just discussed. 

7 Discussion 

In this paper we have performed a detailed analysis of chiral ring partition functions for 
D = 4 superconformal field theories (and analogous BPS states in D = 3, 6) . In particular, 
we have shown how the knowledge of the generating function for the finite N partition 
function, aided with a suitable combinatorial interpretation, can be used to explicitly write 
down the finite N partition functions. These capture the physics of the operators generating 
the chiral ring away from the planar limit. We have also presented detailed analysis of the 
density of states and discussed the notion of typical states for these field theories. 

The finite N result derived above should be useful to understand the behaviour of \- 
BPS and bosonic §-BPS operators in M = 4 SYM, whose finite N partition functions are 
generated by the grand-canonical partition sum (2.3) with k = 2 and k = 3 respectively 
[13]. 31 One motivation for this analysis was to investigate whether it is possible to come up 
with a simple auxiliary model to understand the dynamics of these operators. For instance 
for |-BPS operators one can use the fact that the system has a representation in terms of 
free fermions [56], which has been exploited to understand the detailed behaviour of these 
states from a dual supergravity perspective [43]. 

A natural question is whether there is a similar simplification for | and §-BPS opera- 
tors? 32 This has been discussed in [37] where the answer is argued to be in the affirmative. 
Focussing on the chiral ring, [37] postulates that the dynamics of |-BPS states can be 
encoded by an 5*0(6) matrix model of commuting matrices (see [59, 60] for conjectured ex- 
tensions to M = 1 field theories). One check of this proposal is that it should be able to 
reproduce the correct spectrum of the states in question. While it is possible in this matrix 
model to write down the ground state, the exact spectrum of excited states is not amenable 
to analytic computation, owing to the complicated eigenvalue interaction potential arising 
from the measure. 33 It has been previously argued in [61] that this particular matrix model 
cannot have a free fermion representation and they provide explicit analysis of the interac- 

31 Similar statements can be made for M2 and M5 brane world- volume theories [30] and our considerations 

can easily be extended to these cases. 

32 See recent discussions in [57, 58] where it is argued that perhaps even the non-supersymmetric states 

related to black hole geometries might enjoy a free fermion description. 

33 This is computed using a flat metric on the space of matrices and rewriting this in terms of the eigenvalues 
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tion terms in question. Our analysis is complementary and demonstrates from a spectral 
viewpoint that the system is not simply governed by free fermion dynamics. In particular, 
our factorisation result shows that the partition function for the finite N ^-BPS chiral ring 
of M = 4 SYM can be written as 

Z^ BFS (x 1 ,x 2 ) = Z^ BPS (xi) Z^ BPS (x 2 ) P N (x 1 ,x 2 ) (7.1) 

where P^(xi, x 2 ) is a symmetric polynomial of order ^N(N — 1) in each X{ with non- negative 
coefficients. This very clearly illustrates that the |-BPS operators can be thought of as 
consisting of two sets of different \— BPS operators together with some interaction governed 
by the polynomial. It would be interesting to understand the implication of the factorisation 
result and use it to decode more detailed properties of the chiral ring in M = 4 SYM. 

The main mathematical result of the paper is a derivation of the asymptotic density of 
states for the multi-variable partition functions in question. A special case of this includes 
a generalization of Meinardus' theorem to the case where the associated Dirichlet series 
has multiple poles. These generalizations provide a concrete algorithmic way to decipher 
the asymptotic growth of states in the superconformal theories. Consider the leading order 
estimate (5.66) - for integers which scale as N 2 (i.e., corresponding to operators whose 
charges scale as N 2 ) this grows as \ogd(N 2 ) = O (^N 1 ^ 1 ^. In four dimensional field theories 
with holographic duals, operators of conformal dimension O (N 2 ) are generically expected to 
be dual to heavy states such as black holes. However, since < 2 the growth of the density 
of states in the chiral ring is not enough to account for black hole entropy (a fact noticed by 
many authors previously cf., [13]). Not surprisingly, similar results hold for superconformal 
field theories in other dimensions; for example |— BPS states of the M2-brane worldvolume 
theory has a growth of states given by log g?m2(^) — O (n 4//5 ) ; here the states of interest should 
have conformal dimension O (iV 3 / 2 ) so obtain a non-trivial back-reaction in the dual AdS 
background. The density of such states grows at most as 

N 6/5 < jy3/2_ While thig regult 

can 

be inferred directly without recourse to Meinardus' theorem or generalizations thereof, our 
analysis provides a useful characterization of the sub-leading terms. Assuming one were to 
be able to construct explicit "small-black hole" solutions (i.e., gravitational solutions which 
incorporate higher derivative corrections) dual to these operators and reproduce the leading 
order growth of the density of states, one can then analyze the sub-leading corrections using 
our formalism. 

Our analysis also touched upon the issue of typical operators in the supersymmetric 
sectors; the charge vectors for these operators lie close to the limit curves we derived. It 
would be interesting to understand the relation between the limit curves and the class of 
typical smooth solutions in supergravity. 

On a more technical side, there are a number of open problems. For instance, we have 
briefly touched upon the asymptotics of the finite N partition sums. While the technology 

and the off-diagonal elements. Integrating out the off-diagonal elements leads to the desired measure factor. 
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we developed can be used to extract the asymptotics of the grand canonical partition sum, 
it would be useful to have more explicit finite N asymptotic formulae. Further, we have only 
discussed the asymptotics for the mesonic operators in TV = 1 field theories. Generically, 
these theories also have baryonic operators and the partition sums receive contributions from 
non-zero baryon number sectors. Another interesting class of partition sums are the ones 
that occur in the free theory - these are typically expressed as matrix integrals [12]. Like 
the baryonic partition sums and the finite N results, these generically are not of infinite 
product form, thereby requiring new ideas to extract precise asymptotic formulae. It would 
be interesting to develop the technology to determine the asymptotics of these more general 
partition sums. 
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